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The study of nuclide production and its properties in the N=126 neutron-rich region is the frontier and hot 
topic in nuclear physics and astrophysics research. The upcoming High energy FRagment Separator (HFRS) 
at the High Intensity heavy-ion Accelerator Facility (HIAF), an in-flight separator at relativistic energies, is 
characterized by high beam intensity, large ion-optical acceptance, high magnetic rigidity, and high momentum 
resolution power. It provides an opportunity for the study of the production and properties of neutron-rich nuclei 
around N=126. In this paper, an experimental scheme is proposed to produce the neutron-rich nuclei around 
N=126 and simultaneously measure their mass and lifetime based on the HFRS separator, and the feasibility of 
this scheme is evaluated by the simulations. The results show that under the high resolution optical mode many 
new neutron-rich nuclei approaching the r-process abundance peak around A=195 can be produced for the first 
time, and many nuclei with unknown mass and lifetime can be produced with high statistics. Using the time- 
of-flight corrected by the measured dispersive position and the energy loss information, the cocktails produced 
from the SPb fragmentation can be unambiguously identified. Moreover, the masses of some neutron-rich 
nuclei in the vicinity of N=126 can be measured with high precision using the time-of-flight magnetic rigidity 
technique. This indicates that the HIAF-HFRS facility has potential for the production and property research of 
neutron-rich nuclei around N=126, which is of great significance for expanding chart of nuclides, developing 


nuclear theories, and understanding the origin of heavy elements in the universe. 


Keywords: HFRS, fragmentation, neutron-rich nuclei around N=126, mass measurement, lifetime 


I. INTRODUCTION 


The study of nuclide production and its properties in the 
N=126 neutron-rich region is of great significance to expand 
the nuclear landscape, revealing the evolution of shell struc- 
ture, and understand the astrophysical nucleosynthesis. Ac- 
cording to the theory of nuclear astrophysics, approximately 
half of the nuclei in nature heavier than iron are considered to 
be produced by the rapid neutron capture process, also known 
as the r-process [1-3]. The properties of N=126 neutron-rich 
nuclei, such as mass and lifetime, play a crucial role in un- 
derstanding the r-process abundance peak around A=195 [4— 
6]. However, reaching this region experimentally is difficult 
because of the low production cross-sections and the great 
challenge of separation and identification. The lack of exper- 
imental data leads to significant uncertainty in the predicted 
abundance patterns [7-10]. Therefore, the production and 
property research of neutron-rich nuclei around N=126 is of 
exceptional importance. 

In the last decades, the production of N=126 neutron- 
rich nuclei in the laboratory has been a challenging problem. 
Multi-nucleon transfer (MNT) reactions at near-barrier ener- 
gies are expected to be a powerful technique for the synthesis 
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and study of nuclides in the neutron-rich N=126 region [11-— 
14]. Several facilities based on MNT reactions have been con- 
structed or are in planning or under construction around the 
world, such as KISS at RIKEN [15], IGISOL at JYFL [16], 
the N=126 factory at ANL [17], and INCREASE at GSI [18]. 
Some neutron-rich isotopes of Pt, Ir, and Os, which were pro- 
duced as target-like fragments using the 36 e4!°8pt reaction 
system, have been successfully extracted using the KISS fa- 
cility [19, 20]. However, to reduce the plasma density in 
the gas cell induced by the primary beam and elastic prod- 
ucts [15, 21, 22], the experiments using these facilities are 
usually carried out at low primary beam intensities. More- 
over, thin reaction targets must be used due to the low beam 
energy in the experiments. Low beam intensity and thin re- 
action targets are not beneficial for obtaining high fragment 
yields. To increase the available primary beam intensity, a 
new design using a gas-filled solenoid to suppress unwanted 
elastic particles has been proposed in future facilities such as 
NEXT at Groningen [23] and KISS-II at RIKEN [24]. 


Except for the MNT reactions at near-barrier energies, the 
experimental results indicate that the projectile fragmenta- 
tion reaction of 7°8Pb or PSU at relativistic energies may 
be an effective method for producing heavy neutron-rich nu- 
clei in this region [25, 26]. But this method requires con- 
ditions of high beam intensity, relativistic beam energy, and 
high performance of in-flight separator to be met simulta- 
neously. High beam intensity is a necessary conduction for 
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producing low cross-section products. Relativistic energy 
is required to reduce the number of populated ionic charge 
states, thereby increasing the collection efficiency. High per- 
formance separators are crucial for the separation and iden- 
tification of heavy neutron-rich fragments [27]. The future 
Super-FRS in-flight facility at FAIR [28] meets these condi- 
tions. Searching new isotopes and studying their properties 
in the region along the N=126 line below 7°’Pb have been 
proposed as one of the main physical objectives during FAIR 
Phase-0 experiments [29]. 

High energy FRagment Separator (HFRS) [30, 31], an in- 
flight separator at relativistic energies, is now under con- 
struction at the High Intensity heavy-ion Accelerator Facil- 
ity (HIAF) [32, 33] and will start operation a few years from 
now. It is also characterized by large ion-optical accep- 
tance, high resolution power, high magnetic rigidity, and ex- 
cellent particle identification, just like the Super-FRS separa- 
tor. In combination with the HIAF accelerator facility, which 
will provide beams up to 34 Tm corresponding to 7*8U*>* 
ions of 833 MeV /nucleon and with an intensity as high as 
1 x 10!" ions per pulse, the neutron-rich nuclei around N=126 
from projectile fragmentation at relativistic energies can be 
produced and purified by using the HFRS separator. This pro- 
vides a new opportunity for studying the properties of N=126 
neutron-rich nuclei and understanding the third abundance 
peak in the r-process. 

In this paper, we mainly concentrate on the production 
of N=126 neutron-rich nuclei and the feasibility of studying 
their properties using the HIAF-HFRS facility. The measure- 
ment methods for nuclear properties such as mass and life- 
time, as well as the layout of related experimental setups, will 
be introduced first. Then the simulation results will be pre- 
sented to check the feasibility. 


Il. EXPERIMENTAL SETUPS AND METHODS 


Nuclear masses and lifetimes are of fundamental impor- 
tance for r-process simulations. Using the HFRS separator 
at the HIAF facility, some neutron-rich nuclei in the region 
along the N=126 line below 7°°Pb could be produced, and the 
properties of these nuclei, such as mass and lifetime, could be 
experimentally measured. 

The layout of experimental setups of the HFRS separator 
is shown in Fig. 1. A relativistic 28Pb or PSU primary beam 
will be extracted from the HIAF accelerator facility in a slow 
extraction mode and implanted into the carbon production tar- 
get at the entrance of the HFRS separator. The fragments of 
interest produced via the fragmentation reaction will be col- 
lected and separated using the Bp-AE-Bp method [34, 35], 
in which the analysis of the magnetic rigidity (Bp) is com- 
bined with the energy loss in an achromatic degrader (AF) 
at the pre-separator stage. The achromatic degrader will be 
placed at the PF2 dispersive plane. The unreacted primary 
beam will be intercepted in the beam dump systems installed 
after each dipole magnet in the pre-separator according to Bp 
deviation between the primary beam and the desired isotope. 

The main-separator will be used as a spectrometer. At 
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this stage, the cocktail products will be identified by combin- 
ing the AE-TOF-Bp method and the isomer tagging tech- 
nique [36-38]. The energy loss (AF), time of flight (TOF), 
and the Bp are measured and used to determine the atomic 
number Z and the mass-to-charge ratio A/Q of the frag- 
ments. By detecting delayed y-rays emitted from short-lived 
isomeric states of certain fragment, the cocktails can be un- 
ambiguously identified. The schematic layout of the particle 
identification setups can be seen in Fig. 1. A multiple sam- 
pling ionization chambers (MUSIC) placed at the MF6 will 
be used for the AE measurement [39, 40]. The TOF infor- 
mation will be obtained from the TOF detectors installed at 
the PF4 and MF6 foci with a flight path of 118.03 m [31, 41]. 
The Bp value will be determined from the position measure- 
ments at the MF4 dispersive plane. At the final focal plane of 
the HFRS a decay detector array, which is consists of a sili- 
con array stopper and surrounded by an array of high-purity 
germanium detectors, will be used for the isomer tagging of 
the selected fragments. 

In the AE-TOF-Bp method, the mass-to-charge ratio 
A/Q of a fragment is determined using the equations as fol- 
lows: 


Top = a) 
U 

A Bo 

=R, 2 

5 = wa (2) 


Here L is the length of flight path, v is the velocity of the 
fragment, and y is the relativistic Lorentz factor. The frag- 
ment mass m can then be expressed as, 


cTOF 
L 


m=Q P (1, 6) 
where c represents the velocity of the light. From this equa- 
tion, the nuclear mass can be measured from the measured 
TOF and Bp while performing particle identification. This 
nuclear mass measurement method is called the Bp-TOF 
method [42, 43]. Usually, the precise Bp determination is 
achieved by the trajectory reconstruction method [44]. As- 
suming negligible the object spot size, the Bp values can be 
calculated simply from the following formula, 


TD 


Bp = Bpo[l + tala)” 


(4) 
here Bpo represents the central magnetic rigidity. The zp and 
(x|6) denote the horizontal position and the momentum dis- 
persion at the dispersive plane, respectively. The Bp-TOF 
method has the characteristics of simple equipment, high 
measurement accuracy, low requirement for fragment yield, 
and short measurement time, making it particularly suitable 
for measuring the mass of short-lived nuclei with very low 
yield far from the stability line. Using this method, the masses 
of some very neutron-rich nuclei have been accurately mea- 
sured for the first time, such as the masses of 48-49 Ar and 
557§¢ measured with the combination of the A1900 sepa- 
rator and S800 spectrometer at NSCL [45, 46], the masses of 


156 
157 
158 
159 


160 


161 


162 
163 
164 


16 


a 


166 
167 
168 
169 
170 


17 


172 
173 
174 
175 
176 
177 
178 
179 
180 
181 
182 
183 
184 
185 
186 
187 
188 
189 
190 


19 


192 
193 
194 
195 
196 


197 


Pre-Separator 


| 


Production target Slits TOF1 


Main-Separator 


TOF2 AE 


HHE- 


Position Position 


Position detector 


Fig. 1. (Color online) Schematic layout of experimental setups of the HFRS separator. 


3353 1Ca, a de. 60=62T, and -64y measured with the com- 
bination of the BigRIPS separator and the SHARAQ spec- 
trometer at RIKEN [47, 48]. 

From Eqs. (3) and (4), the mass resolution om /m of the 
Bp-TOF method can be deduced as, 


(Z2)? _ (2 2 (Veror 2 
m Bp TOF (5) 
2 
S Orp 2 Y OTOF \2 
=le Tor 


where o; is the standard deviation of measured values of 
”i”. To improve the mass measurement accuracy, the com- 
mon method is to simultaneously improve the position reso- 
lution of the position detector placed at the dispersive plane, 
the time resolution of the TOF system, and the momentum 
dispersion. For position measurement at the dispersion fo- 
cal plane, detectors with thin material thickness should be 
used to reduce the impact of beam energy loss and multi- 
ple scattering on the beam trajectory. A position-sensitive 
micro-channel-plate (MCP) tracking detector with a reso- 
lution of ~0.5mm (o) [49] and a low-pressure delay-line 
parallel-plate avalanche counter (DL-PPAC) with a resolution 
of ~0.43 mm (ø) [50] were used for the Bp-TOF mass mea- 
surement experiments at NSCL and RIKEN, respectively. Us- 
ing a typical resolution of ~0.5 mm (ø) of the position detec- 
tor at the dispersive plane, combined with a high-resolution 
optical mode with a momentum dispersion of ~10 cm/% spe- 
cially designed for the HFRS main separator, the uncertainty 
oBp/ Bp can be estimated to be about 5 x 1075 

For the TOF measurement, the lower the ion energy, the 
longer its flight time, which is more conducive to improving 
TOF measurement accuracy. However, for N=126 neutron- 
rich heavy nuclei, low ion energy can cause changes in the 
charge state population passing through the materials placed 
at the foci, thereby affecting collection efficiency. Fig. 2 (a) 
shows the fractions of fully stripped ions in the equilibrium 
charge state after passing through an Al degrader for the el- 
ements in the N=126 region as a function of the kinetic en- 
ergy behind the degrader. It was calculated with the Global 
code [51]. One can observe that the fractions of fully stripped 
ions decrease with increasing charge number and decreasing 
kinetic energy. When the ion energy behind the Al degrader 
is greater than 350 MeV/nucleon, even for the heaviest Hg 
element, the proportion of fully stripped ions is higher than 
50%. This proportion is acceptable for the mass measure- 
ment experiments. The flight path length of the HFRS is 
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118.03m. Fig. 2 (b) shows the effect of the time resolu- 
tion of the TOF system on mass measurement accuracy un- 
der different ion energies. For an ion with a kinetic energy of 
350 MeV /nucleon, the contribution of TOF resolution to the 
final mass resolution is ~9.9 x 10~° under a system resolu- 
tion of ~30 ps (a). This TOF resolution requirement can be 
achieved experimentally. For example, at NSCL, two plastic 
scintillator detectors readout by photomultiplier tubes were 
used for the TOF measurement, and the resolution was mea- 
sured to be ~30 ps (o) with primary beam tests [52]. An up- 
graded TOF system for Bp-TOF mass measurement exper- 
iments has been developed at the NSCL and achieved a better 
time resolution (a) of 7.5 ps [53]. In addition, two CVD dia- 
mond detectors formed a TOF system used in RIKEN mass 
measurement experiments, with a system resolution of 27 ps 
(a) [54]. 


Based on the contributions of position measurement ac- 
curacy (~9.9 x 1075) and TOF measurement accuracy 
(~5 x 1075) obtained above, a mass resolution of ~1.1 x 
1074 can be obtained from the Eq. (5), which corresponds 
to Om ~460 keV for ~2000 statistics events and neutron-rich 
nuclei with mass number ~200 near N=126. This accuracy 
of the mass excess is generally sufficient to reveal the evolu- 
tion of shell structure and constrain mass models far from the 
stability [55]. 


Additionally, the mass shifts contributed by isomers will 
be estimated from the isomer measurement using the de- 
cay detector array at the final focal plane of the HFRS. The 
fragments will be stopped in the silicon array stopper after 
the mass measurement and the germanium detectors installed 
close to the stopper will be used to detect the isomeric y- 
rays. From the measurement of the time elapsed between 
the implantation and the subsequent decay, the isomeric life- 
time of the implanted ion can be obtained by correlating the 
particle identification of the HFRS. Simultaneously, using 
the stopper composed of multiple highly-pixelated Double- 
Sided Silicon Strip Detectors (DSSDs) such as AIDA [56] or 
WAS3ABi [57], the implanted nuclei and decay emitting 8- 
rays can be directly correlated within each pixel of the detec- 
tor, providing a direct measurement of the -decay lifetimes. 
This is critical to our understanding of the third abundance 
peak in the r-process. 
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Fig. 2. (Color online) (a) Fractions of fully stripped ions in the equilibrium charge state after passing through an Al degrader for the elements 
in the N=126 region calculated with the Global code. (b) The effect of the time resolution of the TOF system on mass measurement accuracy 


under different ion energies. 


Ill. SIMULATION AND ANALYSIS 


With the HFRS separator at the HIAF facility, many 
neutron-rich nuclei around N=126 could be produced, and 
their properties, such as mass and lifetime, could be exper- 
imentally measured. This is of great significance for expand- 
ing nuclide maps, developing nuclear theories, and under- 
standing the origin of heavy elements in the universe. In 
this section, we will first introduce a high resolution opti- 
cal mode of the HFRS specifically developed for Bp-TOF 
mass measurement experiments. Then the production, sep- 
aration, identification, and nuclear mass measurement accu- 
racy will be studied with the Monte Carlo simulation program 
MOCADI [58] using the high resolution ion-optics for the 
neutron-rich nuclei around N=126. 


A. High resolution ion-optics 


For Bp-TOF mass measurement experiments, the pre- 
separator of the HFRS will be operated as a separator, and the 
main-separator will be used as a spectrometer. A high reso- 
lution ion-optics of the spectrometer is necessary in order to 
accurately measure ion magnetic rigidities. Using the param- 
eterized magnetic field distributions [31], the high resolution 
ion-optics is designed with the codes of Winagile [59] and GI- 
COSY [60]. Fig. 3 shows the beam envelops with the initial 
beam spot sizes of X=+1 mm and Y=+1.5mm. To improve 
the momentum resolution, compared to normal mode [31], 
the horizontal angular acceptance is decreased to +5 mrad 
and the momentum acceptance is reduced to +0.2%. The 
maximum magnetic rigidity is also decreased from 25 Tm to 
15Tm. The vertical angular acceptance remains unchanged 
at 25 mrad. The momentum resolving power at the PF2 and 
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MF4 is 1270 and 7440, respectively, for an emittance of 57 
mm mrad and a horizontal beam spot size of +1 mm. The 
momentum dispersion at the MF4 is 12 cm/%. With this dis- 
persion, the uncertainty of the 7g,/Bp can be estimated to 
be about ~4.17 x 10~° using a typical resolution of ~0.5 mm 
(a) of the position detector at the dispersive plane. 


B. Production and separation 


To produce neutron-rich nuclei around N=126 by fragmen- 
tation reaction, the available projectiles mainly include $Pb 
and $U. In this simulation, both projectiles will be used. 
For a specific nucleus of interest, whether to use Pb fragmen- 
tation or U fragmentation will be determined based on the 
calculated yields. The Booster Ring (BRing) synchrotron of 
the HIAF facility has a maximum magnetic rigidity of 34 Tm 
and can accelerate ?°°Pb*!* and *3°U*>* ions up to 850.74 and 
833.15 MeV/nucleon, respectively. The beam intensities of 
lead and uranium are as high as 1.1 x 1011 and 1.0 x 101" ions 
per pulse, respectively. The typical beam extraction time in 
the slow extraction mode is 3s with a repetition period of 
13s. The corresponding beam spot in the X and Y directions 
on the production target are assumed to be of Gaussian dis- 
tributions with standard deviations of 0.4mm and 0.6mm, 
respectively. 

Taking the production of Au ions as an example, the 
production and separation ability for the neutron-rich nu- 
clei around N=126 will be simulated and studied using the 
high resolution ion-optics of the HFRS. In the simulation, a 
graphite target installed at the PFO focal plane will be used 
as the production target. The thickness of the graphite tar- 
get will be set to 4.4 g/cm? for both Pb and U fragmentation 
reactions, which corresponds to 50.1% of the Pb range and 
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Fig. 3. (Color online) Beam envelopes for high resolution ion-optics. The blue line and the black line correspond to an emittance in the 
X-direction and Y-direction of 57 and 37.57 mm mrad, with an object size of +1 mm and +1.5 mm at the PFO, respectively, while the 


dispersion line (green) represents a momentum deviation of +0.2 %. 


56.5 % of the U range, respectively. In addition, a 75 mg/cm? 
niobium foil is placed behind the production target to achieve 
efficient electron stripping in both reactions. More than 75 % 
of the "Au ions will be fully stripped, which is important for 
high transmission. 


To purify the fragments of interest, an achromatic Al de- 
grader placed at the PF2 dispersive plane will be used. As 
mentioned above, when the ion energy behind the Al de- 
grader is greater than 350 MeV /nucleon, the fraction of fully 
stripped ions is higher than 50% for the neutron-rich nuclei 
around N=126. So the thickness of the degrader will be opti- 
mized to make the energy of fragments of interest greater than 
350 MeV /nucleon. The center thickness of the Al degrader 
is set to 44% and 36.8% of the 7Au range for the Pb and 
U fragmentation reactions, respectively. This thickness com- 
bination of the production target and the degrader will ensure 
that the 7Au ion has approximately 356 MeV /nucleon of 
kinetic energy after passing through the degrader. The Global 
calculations show that the proportion of fully stripped 7Au 
ions at this energy is as high as 57%. Considering the pro- 
portion of fully stripped ions, the HFRS magnetic rigidity is 
set according to the fully stripped ones in both reactions. 


Fig. 4 shows the simulated transmissions of the 7“Au frag- 
ment from the Pb and U fragmentation reactions as a function 
of the separator length. A significant decrease is observed 
near the PF1 and PF2. This is mainly caused by the momen- 
tum deviation. At the PF2 focal plane, the momentum slits 
are set to allow fragments with +0.2 % momentum deviation 
to pass through. In addition, the transmission of the 204 Au 
fragments produced by U fragmentation has a greater reduc- 
tion due to the larger momentum deviation compared to the 
ones of Pb fragmentation. This is because in the U fragmen- 
tation reaction more nucleons need to be wiped out to pro- 
duce %*Au, which results in a greater momentum deviation of 
204 Au fragments according to the Goldhaber model [61]. The 
opening of the mass slits placed at the PF4 focal plane can 
be adjusted based on the desired transmission and the number 
of nuclei required for mass measurements. In the simulation, 
the mass slits are set to +5mm. They also cause a decrease 
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in the transmission due to the charge state population of frag- 
ments. At the final focal plane, the transmissions of 204 Au 
ions are 7 % and 2% for Pb fragmentation and U fragmenta- 
tion, respectively. These transmissions are significantly lower 
than the ones under the normal optical modes described in 
Ref. [31], mainly due to the smaller momentum and horizon- 
tal angular acceptances of the high resolution optical mode. 


The yields and purities of 204 Au ions produced from the Pb 
and U fragmentation reaction are estimated and demonstrated 
in Fig. 5 (a) and (b), respectively. The area of the isotopes in 
the N — Z plane represents the corresponding yields at the 
MF6 focal plane. The yields are obtained from the product of 
the production cross section, the number of target nucleons 
per unit area, and the transmission. The parametrized for- 
mula FRACS can provide a good description for the produc- 
tion cross section estimations [62]. For the number of target 
nucleons per unit area, it is calculated based on the thickness 
of the production target. The transmission is obtained from 
the MOCADI simulations. In addition, to assess the purity, 
nuclei located in the region of Z + 10 and N + 10 around 
24 Au nucleus are selected in the simulation. The purity is 
defined as the ratio of the yield of the fragment of interest to 
the total. From Fig. 5, one can observe that for 204 Ay ions, the 
Pb fragmentation reaction has a yield of up to 9.09 x 10* ppp 
and a purity of about 3.69%. The yield and purity of °” Au 
ions from the Pb fragmentation are higher than those of the U 
fragmentation reaction. Therefore, in future experiments, the 
Pb fragmentation reaction will be selected to produce A 
ions. 


Using the calculated cross sections and the simulated trans- 
missions, the yields of the neutron-rich nuclei around N=126 
are estimated and shown in Fig. 6. In the estimations, the 
primary beam, the production target, and the settings of the 
HFRS are similar to the examples of production and purifica- 
tion for “Au ions. The production cross sections are calcu- 
lated with the FRACS formula. Based on the aforementioned 
simulation results, the transmissions of fragments from Pb 
and U fragmentation reactions are fixed at 7% and 2%, re- 
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Fig. 4. (Color online) Simulated transmissions of the fragmentation product 7“Au as a function of the length of the HFRS separator under 


the high resolution mode. 
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Fig. 5. (Color online) Purification quality of the HFRS for the “Au ions produced by (a) Pb fragmentation and (b) ***U fragmentation 


under high resolution optical mode. 


spectively. The figure presents the results of a 5 —day experi- 
ment, showcasing the maximum yields of fragments of inter- 
est produced through the fragmentation reactions of Pb or U. 
It is obvious that many new neutron-rich nuclei approaching 
the r-process can be produced for the first time, alongside a 
substantial yield of neutron-rich nuclei with unknown mass 
and lifetime. This indicates that the HIAF-HFRS facility has 
great potential for the research on the production and proper- 
ties of neutron-rich nuclei around N=126. 


C. Identification and mass measurement 


Taking the production of “Au ions by Pb fragmentation 
as an example, the particle identification and the accuracy of 
mass measurement using the Bp-TOF method will be stud- 
ied in this section. 

As mentioned above, the particle identification will be 
achieved by the AE-TOF-Bp method. In the simulation, the 
influence of material thickness of TOF and position detec- 
tors on ion transmission and identification is ignored. This is 
because they have a thinner thickness compared to the range 
of the fragments of interest. The thickness of the MUSIC 
energy loss detector is equivalent to 300 wm thick Silicon. 
The time resolution (o) of the TOF system is assumed to be 
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30 ps. The position resolution (ø) of the position detector at 
the dispersive plane and the energy resolution (øo) of the AF 
detector are set to 0.5 mm and 0.4 %, respectively. To obtain 
unambiguous particle identification, it is essential to correct 
the TOF values with the measured dispersive position infor- 
mation at the MF4 focal plane. This can remove the magnetic 
rigidity dependence from the measured TOF spectra. A two 
dimensional correlation between the TOF and the horizon- 
tal position x at the MF4 dispersion plane for Au isotopes is 
shown in Fig. 7 (a). For all nuclides, the deviation between 
their TOF and the central time-of-flight TO Fo can be cor- 
rected using a linear function of the dispersive position x, 

TOF) = TOF — ka. (6) 

Here & represents the slope of the fitted linear function. 
The corrected TOF vs. x spectrum is shown in Fig. 7 (b). 
The original and corrected TOF spectra for Au isotopes are 
compared in Fig. 7 (c). It is obvious that the corrected TOF 
has better resolution, which is beneficial for particle identifi- 
cation. 

The two dimensional correlation spectrum of the measured 
AE and the corrected TOF can provide unambiguous parti- 
cle identification, as shown in Fig. 8. By using the isomer tag- 
ging method, the cocktails can be unambiguously identified. 
In this experimental setting, the charge states between fully 
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Fig. 6. (Color online) Estimated yields of the neutron-rich nuclei around N=126 using the FRACS formula and the simulated transmissions. 
The maximum fragment yields produced by the fragmentation of Pb or U with a 5 —day beam time are shown. The mass and lifetime 
information is extracted from NUBASE2020 [63]. 
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Fig. 7. (Color online) (a) Two dimensional correlation between the TOF and the horizontal position x of the dispersion plane for Au isotopes. 
(b) Corrected TOF vs. x spectrum. (c) Corrected TOF spectra (filled histograms) compared with raw ones (histograms with red lines). 


431 Stripped ions and helium-like ions are observed. As the frag- 432 ments of interest are around A = 2.58 x Z, the hydrogen-like 
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ions 4Z'4~* general appear between the fully stripped ions 
with “°Z7* and ^Z, and the helium-like ions “Z'@~?* 
appear between the ions with 4*°Z7* and 4*°Z* in the TOF 
spectra. This rule is helpful for particle identification. In ad- 
dition, the degrader installed at the PF2 dispersive plane will 
change the charge state population of the fragments passing 
through it. The nuclides circled with red dash line remain 
unchanged in the charge state before and after the degrader. 
The nuclides that have been stripped of one electron by the 
degrader fall into the circle with black dash line. Under the 
same magnetic rigidity, these nuclides have higher velocities 
and shorter flight times. On the contrary, the nuclides that 
capture an electron from the degrader have lower velocities 
and longer flight times, as shown in the circle with green dash 
line in Fig. 8. 

Using the obtained particle identification spectra, we will 
study the mass measurement precision of neutron-rich nu- 
clei around N=126 using the Bp-TOF method on the HFRS. 
The data analysis methodology is basically consistent with 
the Refs. [64, 65]. Typically, nuclei with known masses are 
used to calibrate the relationship between the time-of-flight 
and mass-to-charge ratios. This can help to remove a lot of 
uncertainties in the mass measurement experiment. Among 
the identified cocktails, ten fully stripped nuclides with well- 
known masses according to the 2012 Atomic Mass Evalua- 
tion (AME2012) [66] were selected as calibrants. In addition, 
these chosen nuclides only have known isomers with excita- 
tion energies below 454 keV, which corresponds to the typ- 
ical mass accuracy using the Bp-TOF technique. The cali- 
brants and their atomic mass excess values with uncertainties 
are listed in the left columns of Table 1. The right columns of 
Table 1 list the nuclides to be measured and their mass excess 
values from the literature [66]. 


TABLE 1. Nuclides and their atomic mass excess values with uncer- 
tainties used for the calibration and measurement [66]. Isotopes with 
known isomers with excitation energies below 454 keV are shown 


with ™ , and isotopes with estimated atom mass excess are labeled 
with * 
Calibration nuclides Nuclides to be measured 
oe AME2012 om ians AME2012 om 
(keV) (keV) (keV) (keV) 
MO, -32437.2 2.78 = Ir  -25821#* = 1 967* 
5Qs™ -29511.6 60.55 Ir 216117 196% 
Qs -28278.8 10.06 Ir  -19897% 1967 
1m ~~ -99437.9 38.42 Pt -22692.1 25.15 
197m 228265.8 20.12 pr -19627% 1967 
191, -24400.2 41.06 Au  -23143.5 3.08 
9pm 127390.4 2.22 4 Au* = -20650* 200# 
0p, = -26600.9 20.12 Au” -18770% 1967 
pt  -23740.9 50.1 
202 Au  -24353.0 23.29 


The nuclear masses m can be determined from the corre- 
sponding atomic masses ma utilizing the following formula, 


M = Ma — ZMe + Be(Z), (7) 


where Z and me represent the nuclear charge number and 
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electron rest mass, respectively. B.(Z) denotes the total bind- 
ing energy of the extranuclear electrons. It is calculated with 
an approximate formula [67], 


B.(Z) = 14.438Z7-°9 + 1.55468 x 10-°Z° [eV]. (8) 


From the corrected TOF spectra, the Gaussian fitting func- 
tion is employed to extract the centroids and standard devia- 
tions of the TOF's for both calibration and measurement nu- 
clides. Then we can get the mass calibration function f(r, Z) 
according to the relationship between the time-of-flight cen- 
troids 7 and the mass-to-charge ratios m/q with the following 
relation, 
f(t, Z) =m/q = ao + 17 + aot? + a3TZ + a4Z +052’, 

(9) 
where a; are fit parameters. The terms related to Z are neces- 
sary to account for the impact of the energy loss in the wedge 
degrader at the PF2 dispersive plane. The fit parameters in 
Eq. (9) are determined by an iterative x? minimization proce- 
dure, 


ue a are ZF 
(CAME)? + (Ostat)? 


n 

= > (10) 

i=l 
where n is the number of the calibrants, (m/q)i,ame is 
each calibration mass-to-charge ratio from AME2012 [66], 
(came): and (stat); denote the uncertainty of the mass-to- 
charge ratio from the literature [66] and the statistical uncer- 
tainty, respectively. The statistical uncertainty is calculated 
based on the TOF measurement uncertainty, 


2 
Oi TOF 
(Ostat) N, 


a3Z)°, 


(ai 4 (1) 
where g; ror is the standard deviations of the TOF distri- 
bution, and N; is the statistical count for each calibration. 

After obtaining the fit parameters, the masses of these mea- 
surement nuclides can be calculated from the corresponding 
time-of-flight centroids using Eq. (9). The total uncertainties 
in the mass results mainly include statistical uncertainties, fit- 
ting uncertainties, and systematic uncertainties. The statis- 
tical uncertainty can be estimated with Eq. (11). The fitting 
uncertainty o ;, comes from the uncertainty of the calibration 
function parameters. It is calculated from the error propaga- 
tion based on Eq. (9), 
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where oj; is the covariances of the fit parameters. The fit 


can be performed by the TMinuit class of the CERN root 
package [68], which provides an estimation for the covari- 
ance matrix. The systematic uncertainty o,,, mainly comes 
from the velocity difference caused by the change in charge 
state in the wedge degrader and the method employed to cor- 
rect TOF by dispersive position. We can use the method 
of cross validation of the calibration nuclides to evaluate the 
systematic errors. Assuming there are a total of n calibration 
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Fig. 8. (Color online) Two dimensional correlation spectrum of the measured AE and the corrected TOF. Fully stripped fragments, 
hydrogen-like fragments, and helium-like fragments are labeled with black fonts, red fonts, and green fonts, respectively. The black, red, and 
green circles indicate the isotopes stripped of one electron, the isotopes whose charge states unchanged, and the isotopes captured an electron 


at the PF2 degrader, respectively. 


nuclides, we determine the m/gq value for each of the nuclides ss based on the HFRS separator. 


by calibrating the fitting function of Eq. (9) with the remain- 
ing n — 1 nuclides. Then the normalized chi-square value can 
be calculated as, 


ene 1 3 [(m/q)i,ame — (m/q)i, pit]? 
nonm n & (came); + (ostat); + (Opit)? + Sys" 

(13) 
The differences between the masses from the literature [66] 
and values given by the fitting function as a function of the 
mass-to-charge ratio are shown in Fig. 9 for the calibration 
nuclides (circles) and measured nuclides (squares), and the 
atomic mass excess results and the different contributions 
to the total uncertainties for both the calibration and mea- 
sured nuclides are listed in Table 2. For the statistical un- 
certainties, they are estimated using 5 —day statistical counts 
with Eq. (11). The Ir, ?Pt, and 7°Au nuclei in Table 2 
have significant fitting uncertainties. This is because these 
nuclei have larger mass-to-charge ratios than the calibration 
nuclides, and their masses were extrapolated. If more cali- 
bration nuclides can be selected to cover the mass-to-charge 
ratios of the measured nucleus, smaller fitting uncertainties 
can be obtained. The systematic uncertainty is determined 
as 7.73 keV/q to achieve a normalized chi-square value of 
unity in Eq. (13). From Table 2, the total measurement uncer- 
tainties mainly come from the system uncertainties. This is 
mainly caused by the TOF correction. Using more measure- 
ment information, such as position and angle information at 
the PF4 or MF6, may reduce the system uncertainties. The to- 
tal measurement uncertainty of the measured nuclei is better 
than 900 keV. This measurement accuracy can meet the re- 
quirements of some physical research, proving that the prop- 
erties of some neutron-rich nuclei near N=126 can be studied 


TABLE 2. Nuclides and their atomic mass excess values with un- 
certainties obtained from the mass calibration function. Different 
contributions to the total uncertainties are shown. Isotopes with es- 
timated atom mass excess in the literature [66] are labeled with an 


Uncertainties 

Isotope mit (keV) Ctotai (keV) Ges oe on 

(keV) (keV) (keV) 
140s  -33701.6 606.26 106.78 104.37 587.59 
15Qs™ -28961.5 593.94 77.02 39.7 587.59 
Qs -28447.8 595.61 91.99 32.14 587.59 
ém  -29637.7 598.31 26.19 53.73 595.32 
197m  -28296.1 595.85 22.64 10.88 595.32 
19r -24435.8 596.86 41.33 11.57 595.32 
pe™  -27171.2 603.37 12.0 15.49 603.05 
20Pt  -26692.6 603.15 8.94 6.25 603.05 
201Pt  -23709.1 603.23 12.98 6.77 603.05 
202 Au  -24438.1 611.98 7.84 37.54 610.78 
8  .26374.9 596.01 28.11 11.87 595.32 
207  -21921.2 611.06 61.41 123.38 595.32 
201  -18796.2 868.45 118.62 621.08 595.32 
202Pt  -22580.8 608.4 19.88 78.02 603.05 
0p —_-1 8348.8 743.61 35.16 433.65 603.05 
203Ay  -22687.2 613.34 3.65 55.82 610.78 
24 Au  -19643.5 635.88 5.96 176.81 610.78 
05 Au -16541.0 827.99 10.64 558.92 610.78 
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Fig. 9. (Color online) Differences between the fit masses and the values from AME2012 as a function of the mass-to-charge ratio. The circles 
and squares represent the results of the calibration and measured nuclides, respectively. The error bars of the calibration nuclides include 
the literature uncertainty and the total measurement uncertainty. The uncertainty of the measured nuclides is only the total measurement 


uncertainty. 


IV. SUMMARY 


The properties of N=126 neutron-rich nuclei play crucial 
roles in developing nuclear theories and understanding the 
r-process abundance peak around A=195. To produce these 
neutron-rich nuclei while measuring their mass and lifetime, 
an experimental scheme has been proposed based on the 
HIAF-HFRS facility, and the feasibility of this scheme was 
also evaluated by the simulations in this paper. 

For these studies, a high resolution optical mode of the 
HFRS was specially developed. It has a large momentum 
dispersion of 12cm/% at the MF4 dispersive plane, which 
is beneficial for improving the accuracy of magnetic rigid- 
ity measurement in Bp-TOF mass measurement experi- 
ments. Under this high resolution optical mode, the yields of 
the neutron-rich nuclei around N=126 produced from 7°8Pb 
or PSU fragmentation reactions were estimated using the 
FRACS formula and the simulated transmissions. The results 
show that many new neutron-rich nuclei approaching the r- 
process abundance peak around A=195 can be produced for 
the first time, and many nuclei with unknown mass and life- 
time can be produced with high statistics. Moreover, using the 
time-of-flight corrected by the measured dispersive position 
information and the energy loss information, the cocktails 
produced from the 208PbH fragmentation can be unambigu- 
ously identified, and the masses of some neutron-rich nuclei 
around N=126 can be measured with an accuracy better than 
900 keV by using the Bp-TOF technique. Using these new 
mass data, combined with the machine learning based on the 
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Bayesian neural network as described in Refs. [7, 8], the nu- 
clear mass of this region will be further accurately predicted, 
which is crucial for understanding the r-process abundance 
peak around A=195. These simulation results indicate that 
the HIAF-HFRS facility can provide an opportunity for pro- 
duction and property research of neutron-rich nuclei around 
N=126. 

At present, the HIAF-HFRS facility is under construction. 
All devices including magnets, vacuum, power supply, tar- 
gets, degraders, detectors, etc. are expected to be completed 
and tested by the end of 2024 , after which they can be in- 
stalled on-site. Meanwhile, some high performance detectors 
for the TOF and position measurements used in Bp-TOF 
mass measurement experiments have been proposed. For ex- 
ample, the plastic scintillator detector coupled with multiple 
photomultiplier tube readouts and the diamond detector are 
being developed for time measurements, and the position- 
sensitive micro-channel-plate detector with a large active area 
is also being developed for position measurements. These 
conditions may ensure the carry out of future property re- 
search experiments of neutron-rich nuclei around N=126. 
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